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f*"^ | The AdS/CFT correspondence, or more generally the gauge/gravity duality, is a re- 

£S) ■ markable conjecture obtained from superstring theory with various D-brane backgrounds. 

According to this conjecture, a higher-dimensional curved space-time emerges from super- 
symmetric gauge theory in lower-dimensional flat space-time. In the first part of this article, 
we review Monte Carlo studies of U(7V) supersymmetric gauge theories, which confirmed 
the gauge/gravity duality for various observables. In particular, Monte Carlo results for 
thermodynamic quantities enable us to understand the microscopic origin of the black hole 
entropy associated with the dual geometry. We also discuss results for Wilson loops and 
correlation functions, which agree nicely with the predictions from the gravity side. In the 
second part, we review recent developments in a nonperturbative formulation of superstring 
theory, which may be regarded as a counterpart of the lattice gauge theory in QCD. In 
particular, we discuss Monte Carlo results for the Lorentzian matrix model, which suggest 
that (3+l)-dimensional expanding universe emerges dynamically from type IIB superstring 
Oh theory in (9+1) dimensions if one treats the theory nonperturbatively. 
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§1. Introduction 

Monte Carlo calculations based on lattice gauge theory have been playing a 
central role in studying various properties of QCD in a fully nonperturbative manner. 
In this article, we would like to show that similar developments are starting to take 
place in superstring theory, which is a unified theory for all the matters and the 
fundamental interactions including gravity. As a theory of quantum gravity, the 
notion of space-time had to undergo a revolutionary change. It is by now widely 
' appreciated that matrices are the fundamental degrees of freedom of superstring 

theory at the nonperturbative level, and the space-time only emerges effectively as 
a derived concept at low energy or at long distances. 

One of the developments that manifest this idea is the gauge/gravity duality. 
(See Ref. 1) for a comprehensive review.) The original conjecture is known as the 
AdS/CFT correspondence, which was put forward by Maldacena in 1997, 2 ) but it 
was soon generalized to non-conformal cases. 3 ) These dualities are arrived at by 
considering two different descriptions of D-brane backgrounds in superstring theory. 
One is the gauge theory, which describes the open strings with both ends attached 
to the D-branes. The other is a classical solution of supergravity, which describes 
the closed string degrees of freedom in the bulk sourced by the D-branes. Such 
a remarkable statement that relates gauge theory and gravity theory is possible 
precisely because superstring theory naturally contains both of them. 

When one considers the gauge theory at finite temperature, the dual supergravity 
solution can have a geometry with an event horizon, which is characteristic to black 
holes (or "black branes", more generally). It is well known that black holes have 
thermodynamic properties although their microscopic origin has long been a mystery. 
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The gauge/gravity duality provides a very clear and explicit answer to this problem. 
The thermodynamic properties of a black hole can be understood as those of the dual 
gauge theory, which is considered to describe the interior structure of the black hole. 
Monte Carlo calculation of thermodynamic quantities in the gauge theory indeed 
reproduced the black hole thermodynamics. 4 )~ 6 ) In particular, Ref. 5) shows that the 
duality holds including a' corrections, which represent the effects of closed strings 
having finite length. These results demonstrate that the gauge theory describes 
correctly the quantum space-time structure at the center of the black hole. 

The correspondence between gauge theory and gravity was extended to the level 
of operators, which is of particular importance in using the duality to study strongly 
coupled gauge theories by much simpler calculations in supergravity. In particular, 
explicit prescriptions for obtaining Wilson loops 7 -* and correlation functions 8 ) were 
proposed. The predictions obtained by such prescriptions have also been confirmed 
by direct Monte Carlo calculations on the gauge theory side. 9 ) -11 ) 

Monte Carlo studies mentioned above deal with the simplest case of DO-branes, 
which corresponds to one-dimensional U(iV) super Yang-Mills theory (SYM) with 
16 supercharges. In fact, it is possible to extend these works to higher dimensions by 
extending the idea of the large- N reduction 12 ) to a curved space, 13 ) while circumvent- 
ing a well-known problem 14 ) in the original idea. This unconventional regularization 
scheme, as opposed to the lattice, enables calculations respecting supersymmetry 
maximally. Some preliminary results are obtained 15 ) -17 ) for the case of D3-branes, 
which corresponds to four-dimensional N = 4 U(iV) SYM. 

An important aspect of the gauge/gravity duality is that higher-dimensional 
curved space-time emerges from the gauge theory in lower-dimensional flat space- 
time. The extra spatial dimensions are actually described in gauge theory by the 
scalar fields in the adjoint representation of U(iV) gauge group, which are represented 
by N x N matrices. This is an example of "emergent space" , which appears in various 
contexts of string theory. 18 ) One of the important open questions is whether one 
can extend the idea to emergent space-time instead of just space. Here we quote 
a sentence from Seiberg's lecture 18 ) in 2005: Understanding how time emerges will 
undoubtedly shed new light on some of the most important questions in theoretical 
physics including the origin of the Universe. 

The second development we would like to review in this article concerns a non- 
perturbative formulation of superstring theory, which is considered as a counterpart 
of the lattice gauge theory in the case of QCD. In particular, we discuss recent 
Monte Carlo results for a Lorentzian matrix model, 19 ) which actually show that 
the emergent space-time seems to be naturally realized. Back in 1996, Ishibashi, 
Kawai, Kitazawa and Tsuchiya proposed a matrix model, which is called the type 
IIB matrix model, as a nonperturbative definition of type IIB superstring theory 
in ten dimensions. 20 ) The model can be obtained formally from SYM that appears 
in the aforementioned examples of the gauge/gravity duality by dimensionally re- 
ducing them to zero dimension. Thus one obtains ten bosonic matrices and sixteen 
fermionic matrices, which do not have dependence on space-time coordinates. The 
entire space-time is expected to emerge dynamically from the ten bosonic matrices. 

Until quite recently, however, it was common to study the type IIB matrix model 
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after making a "Wick rotation" . This amounts to replacing the Hermitian matrix Aq 
in the temporal direction by Aq = iA±o, and treating the Hermitian matrix A\q on 
equal footing as the matrices Ai (i = 1, ■ ■ ■ 9) in the spatial directions. The Euclidean 
model obtained in this way has manifest SO (10) symmetry, and it is well defined 
as Monte Carlo studies with small matrices demonstrate. 21 -* In fact the partition 
function was proven to be finite for arbitrary matrix size. 22 ) In Ref. 23), perturbative 
expansion around the diagonal configurations A^ = diag(xi M , • • • , xn>) was studied 
and the low-energy effective theory for the diagonal elements was discussed. In 
particular, it was speculated that configurations with the N points {xf, i = 1, ■ ■ ■ , N} 
distributed on a four-dimensional hypersurface in ten-dimensional Euclidean space 
may be favored due to some nontrivial interactions in the low-energy effective theory. 
If that really happens, it implies that the SO(10) symmetry is spontaneously broken 
down to SO (4) and that four-dimensional space-time is generated dynamically. 

Monte Carlo studies of the Euclidean model is difficult due to the sign problem 
since the Pfaffian that appears from integrating out the fermionic matrices is complex 
in general. Monte Carlo studies of the model omitting the phase of the Pfaffian show 
that the SO(10) symmetry is not spontaneously broken. 24 ) In fact the phase of the 
Pfaffian has an effect of favoring lower-dimensional configurations. 25 ) It is expected 
that such an effect can be studied by Monte Carlo simulation in the near future by 
using a new method to overcome the sign problem 26 )' 27 ) 

As an alternative approach to this issue, the Gaussian expansion method was 
proposed. 28 ) Recently, the free energy was calculated by assuming that the SO(d) 
symmetry (2 < d < 7) remains unbroken, and it was found that d = 3 gives the min- 
imum. 29 ) Another important observation from the Gaussian expansion method was 
that the extent of space-time in the extended d directions and that in the shrunken 
(10 — d) directions turn out to have a finite ratio even in the large- A limit. 29 ) While 
these results reveal interesting dynamical properties of the Euclidean model, which 
can also be understood intuitively from the viewpoint of the low-energy effective 
theory, the connection to our real space-time is not very clear. 

Motivated by these results for the Euclidean model, Kim, J.N. and Tsuchiya 19 ) 
studied the type IIB matrix model without making the Wick rotation. The action has 
an SO(9,l) symmetry instead of SO(10). The reason why no one dared to study this 
Lorentzian model before beyond the classical level 30 ) was that the bosonic action 
is not positive definite unlike the Euclidean case, and therefore the system seemed 
to be highly unstable. Moreover, the bosonic action becomes a pure phase in the 
integrand of the partition function as in the path integral formulation of quantum 
field theories in Minkowski space. Therefore it seemed just impossible to get anything 
sensible out of the Lorentzian model without making a Wick rotation. 

On the other hand, it is known from many examples that the Wick rotation is 
subtle in theories including gravity. For instance, the Lorentzian quantum gravity has 
been pursued within the dynamical triangulation approach 31 ) motivated from earlier 
studies of the Euclidean gravity, and the results turned out to be quite different. More 
recently, the worm hole scenario as a solution to the cosmological constant problem 
was reconsidered in the Lorentzian quantum gravity, 32 ) and the results provided a 
consistent picture, which was not available in the original Euclidean formulation. 
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The crucial trick to make the Lorentzian matrix model accessible by Monte Carlo 
simulation is to integrate out the scale factor of the bosonic matrices first,*) which 
essentially converts the phase factor e tSb into a constraint S*b ~ 0. This is possible 
since the action of the type IIB matrix model is homogeneous* *) with respect to 
the matrices. The model one obtains in this way does not have the sign problem 
since the Pfaffian that appears from integrating out fermionic matrices is real in the 
Lorentzian case. 

First of all, Monte Carlo studies confirmed that the Lorentzian matrix model is 
not well defined as it is. It was found that the extents in the temporal and spatial 
directions, which are represented by -^tr (A)) 2 and ^tr (Ai) 2 , respectively, both tend 
to diverge. One therefore has to put cutoffs on these quantities. However, it turned 
out that the two cutoffs can be removed in the large- N limit in such a way that the 
results scale in N. The theory thus obtained turned out to have no parameters other 
than the scale parameter, which can be naturally identified as the string scale. This 
is a highly nontrivial property of the Lorentzian matrix model, which supports its 
validity as a nonperturbative formulation of superstring theory. 

Another important observation in the Monte Carlo studies is that the eigenvalue 
distribution of the matrix Aq representing the time direction extends as one takes 
the large-iV limit explained above. Supersymmetry plays a crucial role here. (For 
the bosonic model, the eigenvalue distribution of the matrix Aq has finite extent 
even in the large- TV limit.) 

Moreover, by making an SU(iV) transformation in such a way that the tempo- 
ral matrix Aq is diagonalized, the nine-dimensional space represented by A^ (i = 
1, ■ ■ ■ ,9) is found to exhibit a sensible "time evolution". In fact, the space remains 
small and SO(9) symmetric from the infinite past until some "critical time", at 
which only three directions start to expand rapidly. This implies that the rotational 
SO (9) symmetry in the spatial directions is spontaneously broken down to SO (3) 
at the critical time, which may be identified as "the birth of our Universe". The 
so-called initial condition problem is not an issue in the present framework since 
even the time evolution is an emergent concept. The mechanism of the spontaneous 
symmetry breaking (SSB) relies crucially on the noncommutativity of space, and it 
seems totally different from the Euclidean case, in which the SSB is considered to 
be caused by the phase of the Pfaffian. In the Lorentzian model, the Pfafnan is real 
as we mentioned above. 

The rest of this article is organized as follows. In section 2 we review the de- 
velopments related to the gauge/gravity duality. In particular, we discuss its direct 
tests based on Monte Carlo studies of super symmetric gauge theories. In section 3 
we review the developments related to a nonperturbative formulation of superstring 
theory. In particular, we discuss how (3+l)-dimensional expanding universe emerges 
from the Lorentzian matrix model. In section 4 we conclude with a summary and 
future prospects. 



*' The same procedure was also used in Ref. 21) for simulating the Euclidean model. 
**' Note that the gauge theories before dimensional reduction do not have this property due to 
the existence of the derivative terms. 
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§2. Monte Carlo studies of the gauge/gravity duality 

The gauge/gravity duality 1 ) is conjectured by considering D-brane backgrounds 
in superstring theory. D-branes are known to be consistent backgrounds in super- 
string theory, and they may be considered as a counterpart of solitons in field theory. 
Note that solitons do not appear in perturbative expansion around the trivial vac- 
uum, and they are considered as nonperturbative objects. Prom that point of view, 
D-branes are expected to capture some nonperturbative aspects of string theory. 
Indeed D-branes played a crucial role in finding the duality web of superstring/M 
theory and in constructing a nonperturbative formulation of superstring/M theory. 

D-brane can extend in p dimensions, 
and it is characterized as a hypersur- 
face on which strings can end on. ("D" 
stands for the Dirichlet boundary con- 
dition imposed at the boundary of the 
worldsheet of an open string.) Let us 
consider an open string attached to the 
D-brane propagating along it. In Fig. 1 
on the left, we describe such a process 
diagrammatically. If one slices the di- 
agram in the orthogonal direction, one 
notices that the same process can be 
viewed as emission of a closed string. 
This is an example of the so-called open- 
string/closed-string duality. Note here 
that an open string and a closed string 
include a gauge particle and a gravi- 
ton, respectively, as massless modes. This observation lies at the heart of the 
gauge/gravity duality. 

In order to formulate the gauge/gravity duality, we need to consider N D- 
branes lying on top of each other and take a particular low-energy limit so that 
the open strings attached to the D-branes and the closed strings in the bulk are 
decoupled. Then one has two independent descriptions of D-branes. One is the 
(p + l)-dimensional U(7V) SYM, which describes the open strings attached to the 
D-branes. The other is a solution to supergravity, which describes the closed string 
degrees of freedom in the bulk sourced by the D-branes. Note that the bulk ten- 
dimensional space-time is curved since D-branes emit gravitons. In order for the 
supergravity description to be valid, one needs to take certain limits on the gauge 
theory side. 

• The so-called 't Hooft large-iV limit with fixed A = <7ym-^- (The string loop 
corrections are suppressed by 1/N.) 

• The large-A limit. (The a' corrections, which are due to strings having finite 
extent, are suppressed by some powers of 1/A.) 

A stronger version of the conjecture claims that the gauge/gravity duality holds 
including the string loop corrections and the a' corrections. If that is the case, one 




Fig. 1. On the left, an open string attached to 
the D-brane is propagating along it. On 
the right, the same process is viewed as 
emission of a closed string from the D- 
brane. 
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can say that the gauge theory, which is well defined for all N and A, actually defines 
superstring theory on the particular curved background nonperturbatively. 

In addition to this perspective, there are various reasons why the gauge/gravity 
duality are considered very interesting. First of all, it is a realization of an old idea 
by 't Hooft, 33 ) which states that the large-iV gauge theory is equivalent to some 
classical string theory, although in those days people may not have anticipated that 
the string theory actually lives in a curved space-time. It is interesting that the 
curved space-time emerges from a gauge theory in a flat space. This aspect of the 
duality is often referred to as the emergent space. 18 ) In the gauge/gravity duality, 
one typically obtains the anti-de Sitter space. If one considers the gauge theory at 
finite temperature, one can have a black-hole-like geometry. 3 )' 34 ) Therefore, one may 
explain the microscopic origin of the black hole thermodynamics in terms of gauge 
theory. One can also use the duality in the opposite direction, and study strongly 
coupled gauge theories, which are relevant to hadron and condensed matter physics, 
from a curved space-time. 

Since the gauge/gravity duality is a strong- weak duality, it is important to study 
gauge theories in the strongly coupled regime. Monte Carlo simulation can be a pow- 
erful tool for such purposes. However, the problem is that the gauge theories we are 
interested in have supersymmetry, which is broken by the lattice. This can be seen 
immediately if one recalls the supersymmetry algebra {Q, Q} oc P^, where the gener- 
ators for translation appear on the right hand side. Since the translational symmetry 
is broken by the lattice regularization, one necessarily breaks supersymmetry. The 
best one can do is to restore supersymmetry in the continuum limit by fine-tuning 
some parameters in the action, which requires a lot of efforts, however. 

Recently there are considerable developments in "lattice supersymmetry" , which 
can be categorized into two classes. One is the construction of lattice actions with 
various symmetries. For instance, one can preserve one supercharge by using the 
so-called topological twist. The other one, which we discuss in what follows, is to 
use a regularization different from the lattice. 

In the case of DO-branes, which corresponds to the supersymmetric gauge theory 
in 1 dimension with 16 supercharges, one can regularize the theory using momen- 
tum cutoff after fixing the gauge appropriately. 35 ) Black hole dynamics have been 
reproduced including the a' corrections, 4 )' 5 ) and the gauge/gravity duality for the 
Wilson loops 9 ) and the correlation functions 10 )' 11 ) has been confirmed. 

Notably, one can extend this approach to 3d and 4d gauge theories by using 
the idea of large- N reduction. 13 ) In the 4d case, the gauge theory becomes su- 
perconformal and the number of supersymmetries enhances from 16 to 32. This 
superconformal theory is interesting on its own right, but it is also studied inten- 
sively in the context of the AdS/CFT correspondence, which is a typical case of the 
gauge/gravity duality. 2 ) The non-lattice simulation of the 4d superconformal theory 
requires no fine-tuning, unlike the previous proposals based on the lattice regulariza- 
tion. 36 ) As we will see, preliminary results for the Wilson loops and the correlation 
functions are promising. 15 )~ 17 ) 
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2.1. Non-lattice simulation of Id SYM with 16 supercharges 

Let us start with the DO-brane case, which corresponds to the Id U(iV) SYM 
with 16 supercharges. The action is given by S = St, + S{ , where 

Sh = 7 L dttT {\{ DX ^)) 2 - \^)^ X M 2 ) , (2-i) 

Sf = 7 Jo ^ ^ {jf ama ~ \*c*{li)*p[Xi,V^ . (2-2) 

The covariant derivative is denoted as D = <9 t — i[A(t), ■ ]. Xj(t) (j = 1, ■ ■ ■ , 9) 
and •Pa(t) (a = 1, ■ ■ ■ , 16) are N x TV Hermitian matrices, and the theory has SO(9) 
symmetry. When we are interested in finite temperature, we impose periodic bound- 
ary conditions on Xj(t) and anti-periodic boundary conditions on \P a (t). Then the 
temperature is given by T = where /? is the extent in the Euclidean time (t) 
direction. The 't Hooft coupling constant is defined by A = g 2 N, which has the 
dimension of mass cubed. The physics of the system is determined only by the di- 
mensionless coupling constant \ e g = Therefore one can take A = 1 without loss 
of generality. With this convention, the low T regime corresponds to the strongly 
coupled regime, which is expected to have the dual gravity description, 3 ) whereas 
the high T regime is essentially weakly coupled, and the high temperature expansion 
(HTE) is applicable. 37 ) 

In non-lattice simulation, 35 ) we introduce an upper bound on the Fourier mode 
as Xi(t) = Yln=-A Xi :n e twnt , where ui = and similarly for the fermions. This 
idea does not work usually because it breaks gauge invariance. (Recall that the 
Fourier mode is not a gauge invariant concept.) However, in Id, one can fix the 
gauge nonperturbatively in the following way. We first take the static diagonal 
gauge A(t) = -gdiag(ai, • • • , ccjv), in which the gauge field is constant in time and 
diagonal. By following the usual Faddeev-Popov procedure, one obtains 



. a a - a>b 
sin 



(2-3) 



a<b 

as a term to be added to the action. The above gauge choice does not fix the gauge 
symmetry completely, and there is a residual symmetry given by 

a a M- a a + 27ru a , -> Xf n _ Va+Vb , K b n -> Vg n - Va+Vb , (2-4) 

which represents a topologically nontrivial gauge transformation corresponding to 
the gauge function g(t) = diag(e*^ 1 *, • • • ,e tuJUNt ). This residual gauge symmetry 
can be fixed by imposing —it < a a < it. One can then introduce the Fourier mode 
cutoff A. Since there is no UV divergence in this Id model, one can take the A — > oo 
limit naively, and one retrieves the original gauge theory with 16 supercharges. 

The system with finite A can be simulated efficiently *) by using the standard 
RHMC algorithm. 38 ^ In particular, the Fourier acceleration 39 ) can be implemented 

*-* Strictly speaking, the Pfaffian one obtains from integrating out the fermions is complex in 
general. However, there is numerical evidence 11 ' that the phase can be omitted without altering the 
results. Complete justification is left for future investigations. 
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without extra cost since we are dealing with the Fourier modes directly as the fun- 
damental degrees of freedom. This is crucial in reducing the critical slowing down 
at large A. (The same theory is also studied using the standard lattice approach. 6 ) 
However, from the results obtained so far, the non-lattice simulations seem to be far 
more efficient in obtaining the continuum limit.) 

Let us first discuss the phase structure that appears when one changes the tem- 
perature. As is well known, the Polyakov line serves as an order parameter for the 
spontaneous breaking of the center symmetry. Fig. 2 (Left) shows the results. 4 ) At 
high temperature the data agree nicely with the HTE 37 ) including the next-leading 
order. As the temperature decreases below T ~ 3, the data start to deviate, and at 
low temperature below T ~ 0.9, the data can be fitted to the characteristic behavior 
of the "deconfined phase" 

(\P\)=exp(-± + b) . (2-5) 

In the temperature regime investigated, we find no phase transition. This is in 
sharp contrast to the bosonic model, 40 ) -42 ) which undergoes a phase transition to 
the "confined phase" at T ~ 0.9. The absence of the phase transition is consistent 
with analyses on the gravity side. 34 )' 43 ) 
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Fig. 2. (Left) The Polyakov line is plotted against T. The dashed line represents the result of 
HTE up to the next leading order for N = 12. 37 ' The dotted line represents a fit to eq. 
(2-5) with a — 0.15 and 6 = 0.072. (Right) The energy (normalized by iV 2 ) is plotted against 
T. The dashed line represents the result obtained by HTE up to the next leading order for 
N = 12. 37 ' The solid line represents the asymptotic power-law behavior at small T predicted 
by the gauge/gravity duality. 



2.2. Black hole thermodynamics from Id SYM 

Let us turn to a quantitative prediction from the gauge/gravity duality. Given 
the dual geometry, one can use Hawking's theory of black hole thermodynamics to 
obtain various thermodynamic relations such as 44 ) 
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Fig. 3. (Left) The deviation of the internal energy j^E from the leading term 7.41 T~ is plotted 
against the temperature in the log-log scale for A = 1. The solid line represents a fit to a straight 
line with the slope 4.6 predicted from the a' corrections on the gravity side. (Right) The internal 
energy j^jE is plotted against T for A = 1. The solid line represents the leading asymptotic 
behavior at small T predicted by the gauge/gravity duality. The dashed line represents a fit to 
the behavior (2-7) including the subleading term with C = 5.58. 



The gauge/gravity duality predicts that this should be reproduced by Id SYM in 
the large- iV limit at low T. 3 ) The importance of this prediction is that, if it is true, 
it explains the microscopic origin of the black hole thermodynamics, meaning that 
the Id SYM provides the quantum description of the states inside the black hole. 
In Fig. 2 (Right) we plot the internal energy, 4 ) which is defined by E = ■^(fiJ 7 ) 

in terms of the free energy T . At T > 3 the data agree with the HTE. 37 ) As one 
goes to lower temperature, the data points approach the solid line, which represents 
the result (2-6) obtained from the lOd black hole. (See Refs. 45) for earlier studies 
based on the Gaussian approximation.) 

From Fig. 2 (Right) alone, it is not clear whether the gauge theory results 
continue to follow the line predicted from gravity at lower T. In fact, simulations at 
lower T are difficult, since one has to increase A proportionally to 1/T, and at the 
same time one has to increase N to avoid the run-away behavior due to finite iV. 4 ) 
Instead of lowering T, we determined the power of the subleading term as 5 ) 



1 / E \ ( T 



14/5 / rp \ 23/5 

C 7777 , (2-7) 



N 2 VAV3 J ^ A l/3 J ^ A l/3 

from gravity. This was derived by considering higher derivative corrections in the 
supergravity action due to the effects of strings having finite extent (a' corrections). 
The coefficient C of the subleading term is calculable in principle, but it requires the 
full information of the higher derivative corrections, which are yet to be determined. 
By using (2-7), however, we can already make a nontrivial test of the gauge/gravity 
duality. 5 ) In Fig. 3 (Left) we plot the discrepancy 7.41T 14 / 5 - E/N 2 against T 
in the log-log scale, which reveals that the power of the subleading term is indeed 
consistent with the predicted value 23/5 = 4.6. In Fig. 3 (Right) we find that the 
data at T < 0.7 can be nicely fitted to the form (2-7) with C = 5.58. Note also that 
the A = 6 data seem to suffer from some finite A effects at low T. From this point of 
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view, we consider that the A = 4 data points at low T in Fig. 2 (Right), which seem 
to be on the curve of the leading order result from gravity, also suffer from finite 
A effects. Now we know that actually the subleading term in (2-7) should be taken 
into account for precise agreement. 

The fact that the internal energy goes to zero as a power of T towards T = is 
closely related to the existence of the threshold bound state in the Id gauge theory. 46 ) 
The leading power-law behavior T 14 / 5 can actually be understood by considering that 
excitations around the threshold bound state have energy of the order of iV _5//9 as 
suggested from the effective Hamiltonian. In the case of Id SYM with four and eight 
supercharges, it is considered that the threshold bound state does not exist. Indeed, 
Monte Carlo studies show that the internal energy decreases faster as ~ e~ cT as T 
goes to zero. 47 ) 



2.3. Schwarzs child radius from Wilson loop in Id SYM 

As another prediction from the gauge / gravity duality, let us consider the Wilson 
loop, which winds around the temporal direction once. Unlike the usual Polyakov 
line, we consider the one involving the adjoint scalar as 



W = — tr V exp 
N 



P 



dt{A(t) + i ni Xi(t)} 



(2-8) 



where n\ is a unit vector in 9d, which can 
be chosen arbitrarily due to the SO (9) 
invariance. This object can be calcu- 
lated on the gravity side by considering 
the minimal surface spanning the loop 
in the dual geometry. 7 ) For the present 
model, the result is given by 9 ) 

-3/5 



2-ko 1 



— Y 

A 1 / 3 / 



where i?Sch is the Schwarzschild radius 
of the dual black hole geometry and 

2/5 

H - — < — — ■ > = 1.89- •• . 



N=4 




N=6 — 




N=8 — b— 




N=14 




N=17 




high T exp. 


















1/T" 

Fig. 4. The plot of (log \W\) for A = 1 against 
y-3/5 The cuto g- ^ ig c h osen as follows: 

A = 12 for N = 4; A = 0.6/T for N = 6, 8; 
A = 4 for N = 14; A = 6 for N = 17. 
The dashed line represents the results of 
the HTE up to the next-leading order for 
N = 14, which are obtained by applying 
the method in Ref. 37). 



(2-10) 

In Fig. 4 we plot the log of the Wilson 
loop 9 ) against T~ 3 / 5 anticipating (2-9). 
Indeed, at low temperature (to the right 

on the figure), we find that the data points can be fitted nicely to a straight line 
with a slope 1.89 in precise agreement with (2-10). The solid line corresponds to 
(log|W|) = 1.89 T -3 / 5 — 4.58, where the existence of the constant term can be 
understood as a' corrections. This result demonstrates that one can extract the 
information of the dual geometry such as the Schwarzschild radius from the gauge 
invariant observable (2-8). 
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2.4. Correlation functions in Id SYM 

One can also predict correlation functions from the gravity side. This was done 
more than ten years ago by Sekino and Yoneya 48 ) applying the Gubser-Klebanov- 
Polykov-Witten prescription 8 ) to the present Id SYM case. 

For instance, let us consider an op- 
erator 

O e = ±Stv(F ij X il ...X ie ) , (2-11) 

where = —i[Xi,Xj] and Str repre- 
sents the symmetrized trace treating F^ 
as a single unit. The two-point corre- 
lation function of this operator is pre- 
dicted as 

<0<(t)<W)> ~ ^ , P=J + 1 

(2-12) 

at A" 1 / 3 < \t\ < A-VSjV 10 / 21 . In Fig. 5 
we plot the two-point correlation func- 
tion 11 ) for £ = 1,2,3,4, which agrees 
precisely with the predicted power-law 
behavior.*) What is rather surprising is 
that the agreement is observed even for 
such small N as N = 3. In particular, the power-law behavior seems to extend to 
the far infrared regime, in which the supergravity calculations become invalid. This 
fact may have important implications on M theory interpretation of the same model. 
See Ref. 11) for more details as well as results for other operators. 

In general, the correlation functions of operators which correspond to super- 
gravity modes on the gravity side show power-law behavior as predicted by the 
gauge/gravity duality. Some operators are predicted to show unusual infrared di- 
verging behavior such as |p| -6 / 5 in the momentum space. The gauge/gravity duality 
is confirmed by Monte Carlo calculations even in such cases. The operators corre- 
sponding to stringy excited modes are also studied on the gravity side, 50 ) and they 
are predicted to have correlation functions with an exotic behavior e~ ct3/ ° . Monte 
Carlo results for these operators are indeed consistent with such a behavior. 

2.5. Extension to higher dimensions based on the large-N reduction 

In this section we discuss how one can extend the works in the previous sections 
to higher dimensions. Respecting supersymmetry becomes more nontrivial in higher 
dimensions, and we use the idea of the large- N reduction, which relates gauge theories 
in higher dimensions to those in lower dimensions in the 't Hooft large- N limit. This, 

As a closely related work, Refs. 49) presented a numerical analysis based on the discrete light- 
cone quantization for the (1 + 1) -dimensional case in the large-iV limit. The two-point correlation 
function of the stress-energy tensor has been calculated, and the expected power-law behavior has 
been confirmed. 
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The cutoff parameters are chosen as /3 = 4 
and A — 16. The straight lines represent 
the power-law behavior predicted by the 
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in particular, enables us to perform Monte Carlo studies of the D3-brane case, which 
corresponds to the four-dimensional N = 4 SYM. 

Here we use a novel large- N reduction, which was proposed in Ref. 13) to study 
M = 4 SYM on R x S 3 . The reduced model in this case is obtained by collapsing 
the S 3 to a point, and it is given by the action 



5pw = — — / dt ti 
9pw J 



X -(p t X M ? - \[X M ,X N } 2 + U>D t * - hp7 M [X M ,¥] 



+^(X l ) 2 + ^{X a f + i^e ijk X t X 3 X k + i^ 7 i23^ 



(2-13) 



where the parameter \i is related to the radius of S 3 as R s s = ^, and the covariant 
derivative is defined by Dt = dt — i[A, ■ ], where A(t) as well as Xm{P) and &(t) 
is an N x N Hermitian matrix. The range of indices is given by 1 < M, N < 9, 
1 < k < 3 and 4 < a < 9, respectively. The model has the SU(2|4) symmetry 
with 16 supercharges. 

The one-dimensional super symmetric gauge theory (2-13) is called the plane 
wave matrix model (PWMM) or the BMN matrix model 51 )*). It is nothing but 
the Id SYM discussed in the previous sections, plus some mass deformation, which 
preserves 16 supersymmetries of the undeformed theory. 

The PWMM possesses many discrete vacua representing multi fuzzy spheres, 
which are given explicitly by 

V V 

Xi = /i0(Lf l) 8l /!/ ) with ^n 7 fc 7 = iV, (2-14) 
i=i i=i 

where l}p are the r-dimensional irreducible representation of the SU(2) algebra 
[L^ r \ L^p] = i€ijkL^\ These vacua preserve the SU(2|4) symmetry, and are all 
degenerate. 

In order to retrieve the planar M = 4 SYM on R x S 3 , one has to pick up a 
particular background from (2-14), and consider the theory (2-13) around it. Let us 
consider the vacuum denned by 

ki = k , nj = n + / — — ^— ■ f° r ^ = lr""> I/ > (2T5) 
and take the large- iV limit in such a way that 

k — > qg , ^oo, v — > oo , with Apw = ^ PW fixed . (2T6) 

v n 

Then the resulting theory is claimed 13 -* to be equivalent**) to the planar limit of 

*^ Properties of this model at finite temperature are studied at weak coupling 52 '' 53 ' and at 
strong coupling. 54 ' 

**' See Refs. 55) for earlier studies that led to this proposal. This equivalence was checked at 
finite temperature in the weak coupling regime. 53 ' It has also been extended to general group 
manifolds and coset spaces. 56 ' 
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M = 4 SYM on R x S 3 with the 't Hooft coupling constant given by 

r3 



Asym = 2vr 2 Ap W (i?53) 3 = ™L±2ZZ. . (2-17) 

n ji 6 



In practice, we use (2-14) with (2-15) as the initial configuration and check that no 
transition to other vacua occurs during the simulation. 

The above equivalence may be viewed as an extension of the large- ./V reduc- 
tion, 12 ) which asserts that the large- ./V gauge theories can be studied by dimension- 
ally reduced models. It is known that the original idea for theories compactified on 
a torus can fail due to the instability of the U(1) D symmetric vacuum of the re- 
duced model. 14 ) This problem is avoided in the novel proposal since the PWMM is 
a massive theory and the vacuum preserves the maximal SUSY. This regularization 
respects 16 super symmetries, which is half of the full superconformal symmetry of 
N = 4 SYM on R x S 3 . Since any kind of UV regularization breaks the conformal 
symmetry, this regularization is optimal from the viewpoint of preserving SUSY. 

Since the parameter gp w in the action (2T3) can be scaled out by appropriate 
redefinition of fields and parameters, we take gp W -/V = 1 without loss of generality 
as in Refs. 4),5),9)-ll). In this convention one finds from eq. (2-17) that the small 
(large) ji region in the PWMM corresponds to the strong (weak) coupling region in 
the 4d N = 4 SYM. 

2.6. Wilson loops in 4d M = 4 SYM 

Let us consider the following type of Wilson loop 

W(C) = ^tvVexp £ds (iAfx»{s) + \x»{s)\ xfe a ) , (2-18) 

where x^(s) = dx d ^ and 9 a is a constant which satisfies 6 a 9 a = 1. The fields A* 
and X^ 4 represent the gauge field and the six scalars, respectively, in 4d M = 4 
SYM on R 4 . Due to the particular way in which the scalars appear, one can obtain 
predictions from the gravity side based on the AdS/CFT correspondence as 7 ) 

lim (W(C)\ = e~ s{c) , (2-19) 

AT-+oo,A S ym->-oo \ /SYM 

where S{C) represents the area of the minimal surface spanning the loop C on the 
boundary of the AdS space. 

For the circular Wilson loop IU(C c i rc ), which is a (l/2-)BPS operator, there is 
an exact result on the gauge theory side, which is obtained by summing up planar 
ladder diagrams 57 ) or by using the localization method. 58 ) The result is given by 



lim 

7V-s>oo 



W(C c « c )) sYu = ^h(V^) (2-20) 

for Asym » 1 , (2-21) 



gV / 2AsYM 



(f) 1/2 (2A S YM) 3 / 4 
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where I\ (x) is the modified Bessel function of the first kind. The result is independent 
of the radius of the circle, which is a consequence of the scale invariance of M = 4 
SYM. At strong coupling it agrees with the result obtained from the dual geometry 59 ) 
S(C c i TC ) = — v / 2Asym- This is an explicit example of the AdS/CFT correspondence. 
We use the exact result (2-20) for arbitrary Asym to test our calculation method. 

The Wilson loop in N = 4 SYM can be calculated in PWMM in the following 
way. When we perform the conformal mapping from R 4 to R x S* 3 , the radial and 
angular directions are mapped to the time and ^-directions, respectively. Therefore, 
an arbitrary loop on a plane in R 4 is mapped to a loop on R x S 3 , which can be 
projected to a great circle on S 3 . Such a Wilson loop can be represented in the 
large- N reduced model as 

W rcd (C) = -^trVexp£ds (iAo^- +iXi^{s) + \x^(s)\ XJ a ) , (2-22) 

where e* (x^(s)) is the dreibein on S 3 . 

The expectation value of this oper- 
ator is related to the average of the orig- 
inal Wilson loop as 60 ^ 

(^ (C) )sYM = ( W »*W) ' ^ 

where (• • • ) on the right-hand side de- 
notes the expectation value in the large- 
N reduced model (PWMM). In the case 
of the circular Wilson loop, the rela- 
tion (2-23) was confirmed by reproduc- 
ing the SYM result (2-20) from the re- 
duced model to all orders in perturba- 
tion theory assuming that non-ladder di- 
agrams do not contribute. 61 ^ 

In Fig. 6 we present our preliminary 
results for the circular Wilson loop. 17 -* 
We have performed the A — > oo extrapo- 
lation using A = 6,8,10,12 assuming that 
finite A effects are 0(1/ A). The extent in the time direction is fixed to (3 = 5. The 
parameters describing the background (2T5) are chosen to be (n,u) = (|,2), and 
we performed an extrapolation to k = oo using the data for k = 2, 3, 4, 5 assuming 
that the finite-k effects are 0(l/k 2 ). We also plot the exact result (2-20). Except 
for the data point at \Asym = 4, the agreement with the exact result is promising. 
Note, in particular, that we already start to observe a bent from the weak coupling 
behavior towards the strong coupling behavior. This is remarkable considering the 
rather small matrix size. We consider this as a result of the fact that our formulation 
respects sixteen supersymmetries. 




VAsym 



Fig. 6. The log of the circular Wilson loop 
normalized by \J Asym is plotted against 
V-^sym- The solid line represents the ex- 
act result (2-20). The dashed line repre- 
sents the behavior (2-21) at strong cou- 
pling, whereas the dotted line represents 
the leading perturbative behavior In (W) ~ 
4 Asym. 
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2.7. Correlation functions in J^d M = 4 SYM 

In this section we consider chiral primary operators (CPOs) as simple examples 
of 1/2 BPS operators in 4d N = 4 SYM, and present Monte Carlo results for their 
correlation functions. 16 ) In particular, we find that the two-point and three-point 
functions agree with the free theory results up to overall constant factors even at 
fairly strong coupling. Moreover the ratio of the overall factors agrees with the 
prediction of the AdS/CFT correspondence. *) 




Fig. 7. (Left) The two-point function (trZ 2 (p) trZ+ 2 (— p)j is plotted in the log scale. 
The curves represent the corresponding free theory results multiplied by 0.919, 
0.799, 0.647 for fx = 4.0, 2.0, 1.3, respectively. (Right) The three-point function 
^tr (x4,X$(p)^ tr (^X^Xq^)^ tr (^X^X^— p)^ ^> is plotted in the log scale. The curves repre- 
sent the corresponding free theory results multiplied by 0.850, 0.716, 0.491 for fi = 4.0, 2.0, 1.3, 
respectively. 

Let us consider the CPOs given by 

of(x) = T ai ... aA tr (X*X* • • • X*>)) , (2-24) 

where T ai ... aA is a symmetric traceless tensor and X^ 4 represents the six scalar s in 
4d M = 4 SYM on R 4 . Thanks to the conformal symmetry, the forms of two-point 
and three-point functions of the CPOs are determined as 

0f( Xl )0f(x 2 )) = c A (of(x 1 )of(x 2 )) , 

I \ I free 

R 4 <„ \/nR 4 /„ \snR 4 1„ \\ //nR 4 /„ \snR 4 1 „ \snR 4 



OX(^m 2 (^)OX(x 3 )) = c AlA2A , {OX(xi)02 2 (x 2 )02 3 (x 3 )) bee ,(2-25) 

where c A and c AlA2A:i axe over-all constants depending on Asym m general, and 
(• • • )frcc denotes the results of free theory. The analysis on the gravity side suggests 63 ) 



C-A 1 A 2 A 3 



^/c Al c A2 c As 



C A\A 2 A Z 



JV->-oo,Asym->-oo \f c A 1 CA 2 CA 3 



= 1 for M; . (2-26) 

Ar-*-oo,A S YM->-0 



*-* There are also Monte Carlo studies of the 4d M = 4 SYM based on matrix quantum mechanics 
of 6 bosonic commuting matrices, 62 ' which give results consistent with the AdS/CFT for the three- 
point functions of CPOs. 
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In order to relate the above operators to those in the PWMM, we first perform 
the conformal mapping*) from R A to R x S 3 . Then the M-point functions of the 
CPO Of" 8 " on R x S 3 are related to those in PWMM as 



J ~2^"' J 



da 



(M) 



2vr 2 



O 



RxS 3 



(M) 
3 



(O 



PW 



(ii)-^!N) , 



(2-27) 
(2-28) 



where we have defined 13 ) 0^{t) = T c 



tr (x ai X a2 



X, 



a A 



it))- 
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We calculate the two-point functions (tr Z 2 {t\) tr Z^ 2 (t-iyj , where Z = ^75(^4 + 

2X5), and the three-point functions (tr (^X^X^iti)) tr ^Xs-Xg^)^ tr (^X^X^t^))^ . 

The CPOs we consider here have A = 2, and the AdS/CFT predicts C222 = ci^ 2 , 
which we test by Monte Carlo calculations. 

The parameters describing the 
background (2T5) are chosen as n = 

= 2,k = 2, which corresponds to 
the matrix size N = 6. The values of 
H we use are ji = 4.0,2.0,1.3, which 
correspond to Asym — 0.55, 4.39, 16.0, 
respectively, in the chosen background. 
Thus we cover a wide range of the cou- 
pling constant. The regularization pa- 
rameters in the t-direction are taken as 
(3 = 5.0 and A = 12 for all cases. 

In Fig. 7 (Left) we plot the two- 
point function**) <^tr Z*(p) tr Z^^-p)^. 

We find that the results agree well — 
up to overall constants depending on 
/j, — with the corresponding free the- 
ory results, which are obtained analytically by just switching off the interac- 
tion terms in the reduced model with the same regularization parameters. In 
Fig. 7 (Right) we show similar results for the three-point function defined by 

(tr (jUX 5 {p)) tr (XXe(0)) tr (x^X 4 {-p)) ). 

We can extract the the overall constants corresponding to C2 and C222 in eq. 
(2-25) from Fig. 7. In Fig. 8 we plot the overall constants obtained in this way 
for three values of /x. The data points represent the mean value of the upper and 
lower bounds. We find that our results for various coupling constants lie on the 



0.6 



Fig. 8. The overall constants corresponding to 
C2 and C222 in eq. (2-25) are plotted in the 
log-log scale. The straight line presents 
the relation C222 = c^ 2 predicted by the 
AdS/CFT. 



straight line which represents the prediction C222 



„ 3 / 2 



from the AdS/CFT. Our 



9 i=-t 



dr z + r A dni 



^ The metrics of R 4 and RxS 3 are related as ds 2 R4 . 

The transformation of the CPOs is given by 0% xSd = e^^O% 
-* The Fourier transform of an operator 0(t) is defined as <D(p) = jj dtO(t) e~ %pt 



fixS 3 ' 



where 
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results therefore suggest that the relation (2-26) holds also at intermediate coupling 
constants. 

§3. Nonperturbative formulation of superstring theory 

Superstring theory not only provides a most natural candidate for a consistent 
theory of quantum gravity but also enables unified description of all the interactions 
and the matters. A crucial problem is that we do not yet have a well-established 
nonperturbative formulation, which would be needed in addressing dynamical issues 
such as the determination of space-time dimensionality. *) 

In the 1990s, there was a remarkable progress in understanding the nonpertur- 
bative aspects of superstring theory based on D-branes. Most importantly, it was no- 
ticed that large- N matrices are the appropriate microscopic degrees of freedom which 
are useful in formulating superstring theory in a nonperturbative manner. 20 )' 65 )' 66 ) 
In particular, the type IIB matrix model was proposed as a nonperturbative formu- 
lation of type IIB superstring theory in ten-dimensional space-time. 20 ) It was also 
realized that the five types of superstring theory in ten dimensions are just differ- 
ent descriptions of the same theory. Therefore, it was speculated that the type IIB 
matrix model actually describes the unique underlying theory, although it takes the 
form that has explicit connection to perturbative type IIB superstring theory. 20 )' 67 ) 

In the type IIB matrix model, the space-time is represented dynamically by the 
eigenvalue distribution of ten bosonic NxN traceless Hermitian matrices. 23 ) We first 
discuss the Euclidean model, which is obtained by making a Wick rotation. Then we 
discuss a recent work on the Lorentzian model, which shows that (3+l)-dimensional 
expanding universe emerges dynamically. 19 ) 

3.1. Definition of the Euclidean model 

The type IIB matrix model can be obtained formally by the zero-volume limit 
of D = 10 SU(iV) pure super Yang-Mills theory. The partition function of the 
Euclidean model is given by 

Z = j dAd&e- Sb - S{ , (3-1) 
Sb = -^tv[A^A v ] 2 , (3-2) 
St = ~tr {* a {Cr») aP [A^}) . (3-3) 

Here A^ (// = 1, ■ ■ ■ , 10) are traceless NxN Hermitian matrices, whereas \P a (a = 
1, ■ ■ ■ , 16) are traceless NxN matrices with Grassmannian entries. The parameter 
g can be scaled out by appropriate redefinition of the matrices, and hence it is just 
a scale parameter. We therefore set g 2 N = 1 from now on without loss of generality. 



See Ref. 64) for a well-known scenario based on string-gas cosmology. 
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The integration measure for and \P a is given by 

N 2 -l 10 iV 2 -l 16 

a=l /i=l a=l a=l 

where A a and are the coefficients in the expansion = Yla^ 1 A^T 0, etc. with 
respect to the SU(iV) generators T a normalized as tr (T a T b ) = \5 ab . 

The model has an SO (10) symmetry, under which A^ and \P a transform as a 
vector and a Majorana-Weyl spinor, respectively. The 16 x 16 matrices T^ are the 
gamma matrices after the Weyl projection, and C is the charge conjugation matrix, 
which satisfies (T^) T = CT^tf and C T = C. 

In general, one can obtain super symmetric matrix models by taking the zero- 
volume limit of pure super Yang-Mills theories in D = 3, 4, 6 and 10 dimensions, 
where the D = 10 case corresponds to the type IIB matrix model. The conver- 
gence of the partition function for general D was investigated both numerically 21 ) 
and analytically. 22 ^ The D = 3 model is ill-defined since the partition function is 
divergent. The D = 4 model has a real positive fermion determinant, and Monte 
Carlo simulation suggested the absence of the SSB of rotational symmetry. 68 ) (Sec 
also Refs. 69), 70).) The D = 6 model and the D = 10 model both have a complex 
fermion determinant, whose phase is expected to play a crucial role 25 )~ 27 )' 71 )' 72 ) in 
the SSB of SO(L>). 

In order to discuss the spontaneous symmetry breaking (SSB) of SO(10) in the 
large- TV limit, we consider the "moment of inertia" tensor 23 )' 73 ) 

T lw = ^tT(A li A u ) , (3-5) 

which is a 10 x 10 real symmetric tensor. We denote its eigenvalues as Xj (j = 
1, ■ ■ ■ , 10) with the specific order 

Ai > A 2 > • • • > Aio . (3-6) 

If the SO(10) is not spontaneously broken, the expectation values (Xj) (j = 1, ■ ■ ■ , 10) 
should be all equal in the large- iV limit. Therefore, if we find that they are not equal, 
it implies that the SO(10) symmetry is spontaneously broken. Thus the expectation 
values (Xj) serve as an order parameter of the SSB. 

3.2. The Gaussian expansion method 

Since there are no quadratic terms in the actions (3-2) and (3-3), we cannot 
perform a perturbative expansion in the ordinary sense. Finding the vacuum of 
this model is therefore a problem of solving a strongly coupled system. Here we 
review the results obtained recently by the Gaussian expansion method. 29 ) See also 
Refs. 28), 74) for earlier works. The application of such a method to large- N matrix 
quantum mechanics was advocated by Kabat and Lifschytz, 75 ) and various black 
hole physics of the dual geometry has been discussed. 45 ) 
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The starting point of the Gaussian expansion method is to introduce a Gaussian 
term Sq and to rewrite the action S = Sb + Sf as 

S=(S + S)-S . (3-7) 

Then we can perform a perturbative expansion regarding the first term (So + S) as 
the "classical action" and the second term (—So) as the "one-loop counter term". 
The results at finite order depend, of course, on the choice of the Gaussian term So, 
which contains many free parameters in general. However, it is known in various 
examples that there exists a region of parameters, in which the results obtained at 
finite order are almost constant. Therefore, if we can identify this "plateau region", 
we can make concrete predictions. It should be emphasized that the method enables 
us to obtain genuinely nonperturbative results, although most of the tasks involved 
are nothing more than perturbative calculations. 76 -* There are some cases in which 
one finds more than one plateau regions in the parameter space. In that case, each 
of them is considered to correspond to a local minimum of the effective action, and 
the plateau which gives the smallest free energy corresponds to the true vacuum. 
These statements have been confirmed explicitly in simpler matrix models. 77 ) 

As the Gaussian action for the present model, we consider the most general one 
that preserves the SU(iV) symmetry. Note, in particular, that we have to allow the 
Gaussian action to break the SO(10) symmetry so that we can study the SSB of 
SO(10). Making use of the SO(10) symmetry of the model, we can always bring the 
Gaussian action into the form 

N 10 N 16 

H=l «,/3=l 

where and A a p are arbitrary parameters. The 16 x 16 complex matrix A a p can 
be expanded in term of the gamma matrices as 

10 

H,v,p=l 

using a 3- form m^ up . 

In practice, we truncate the series expansion at some finite order. Then the 
free energy and the observable depend on the free parameters M M and A a /3 in the 
Gaussian action. We search for the values of parameters, at which the free energy 
becomes stationary by solving the "self-consistency equations" 

and estimate F and observables at the solutions. As we increase the order of the 
expansion, the number of solutions increases. If we find that there are many solutions 
close to each other in the parameter space which give similar results for the free 
energy and the observables, we may identify the region as a plateau. In fact we 
restrict the parameter space by imposing the SO(d) symmetry with 2 < d < 7. The 
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Fig. 9. (Left) The free energy density averaged over the "physical solutions" for each d is plotted 
against d. The horizontal line represents the KNS value / = log 8 — § = 1.32944 . . ., and the 
dotted line connecting the data points is drawn to guide the eye. (Right) The extent of space- 
time R 2 and r 2 in the extended and shrunken directions, respectively, are plotted against d 



after taking the average over the "physical solutions" for each d. 
connecting the data points are drawn to guide the eye. 



The solid and dashed lines 



plateau region identified for each d corresponds to a local minimum which breaks 
the SO(10) symmetry spontaneously down to SO(d). By comparing the free energy, 
we can determine which local minimum is actually the true vacuum. 

For each d, we first obtain the free energy up to the third order as a function 
of the free parameters in the Gaussian action. More precisely, we actually calculate 
"the free energy density" defined as 

/= lim /_Z_-(-31og7V)) , where F = -log Z . (3-11) 

By differentiating the free energy density with respect to the free parameters, we 
obtain the self-consistency equations, which we solve numerically by Mathematica. 

In Fig. 9 (Left) we show the free energy density obtained by averaging over the 
"physical solutions" (See Ref. 29) for more detail.) for each d at order 3. We put 
error bars representing the mean square error when there are more than one physical 
solutions. The horizontal line log 8 — | = 1.32944... represents the result obtained 
from the analytic formula of the partition function conjectured by Krauth, Nicolai 
and Staudacher (KNS). 21 ) 

The result decreases monotonically as d decreases from 7 to 3, and it becomes 
much larger for d = 2. Thus, the SO(3) symmetric vacuum gives the smallest free 
energy density. The (^-dependence of the free energy density is quite analogous to 
the one observed in the six-dimensional case. 78 ) There the value of the free energy 
tends to decrease slightly as one goes from order 3 to order 5. Considering such 
artifacts due to truncation, we speculate that the KNS conjecture actually refers to 
the partition function for the SO (10) symmetric vacuum. 

Let us discuss the results for the extent of space-time represented by the eigen- 
values (3-6). For each of the SO(d) symmetric vacua, the d large eigenvalues 
(1 < j < d) are equal due to the imposed SO(d) symmetry, and we denote the corre- 
sponding value as R 2 . The remaining (10 — d) eigenvalues for each solution turn out 
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to be quite close to each other and we denote the mean value as r 2 . In Fig. 9 (Right) 
we plot the result for R 2 and r 2 averaged over all the physical solutions for each d. 
We put error bars representing the mean square error when there are more than one 
physical solutions. We find that r 2 stays almost constant at r 2 = 0.1 ~ 0.15, which 
seems to be universal for all the SO(d) symmetric vacua with 2 < d < 7. On the 
other hand, the results for R 2 are found to be larger for smaller d. It turned out 
that this behavior is consistent with the constant-volume property, 78 ^ which is given 
by R d r 10 - d » £ 10 with I 2 « 0.38, except for d = 2. 

3.3. The mechanism of SSB in the Euclidean model 

In this subsection, we review the arguments in Ref. 25) which show that the 
phase of the Pfaffian favors d(> 3)-dimensional configurations with Xj (j = d + 
1, ■ ■ • > 10) much smaller than the others. This suggests the possibility that the 
SO(10) symmetry is broken down to SO(d) with d > 3. 

Going back to the definition of the model (3-1), let us first integrate over the 
fermionic matrices which yields 



dip- e~ Si = PfM , (3-12) 
where 

Maa,b/3 = ~i fabc(C T^apA^ (3T3) 

is a 16 (TV 2 — 1) x 16 (TV 2 — 1) anti-symmetric matrix, regarding each of (aa) and 
(b/3) as a single index. The real totally-antisymmetric tensor f a b c gives the structure 
constants of SU(iV). In what follows, it proves convenient to work with an explicit 
representation of the gamma matrices given by 

A = i (72 <8> 02 ® (7 2 <g> (72 , r 2 = i <T 2 <8> 02 ® 1 ® °~1 > A = ^ 02 <8> 02 <8> 1 <8> 03 , 
A = i (7 2 (8) (7i (8) 0~2 <8> 1 , r 5 = « (7 2 (8) <T 3 (g) <7 2 <8> 1 , A = « (7 2 <8> 1 <8> CTi (g) (7 2 , 
A = * o- 2 (8) 1 (8) (7 3 (g) <7 2 , r 8 = j(7i<8)l®ll8ll , r 9 = i(73(8)l(8)l(8)l, 

Ao = 1®1®1®1, (3-14) 

for which the charge conjugation matrix C becomes a unit matrix. Note that A is 
pure imaginary except for Ao i n this representation. Therefore, the Pfaffian is real 
for A±o = (or for the Lorentzian case A±o = — iAq). 

Next, when A 3 = A 4 = ■ ■ ■ = A w = 0, one finds that Pf M(A) = 0. 21 ) Note 
first that (Pf.M) 2 = det M = \ detC/| 16 , where U is an (N 2 - 1) x (N 2 - 1) matrix 
defined as U a b = f a bcX c , where X c = A\ +iA\. Since U a bX b = 0, the matrix U has 
a zero-eigenvalue, and therefore det {7 = 0. 

Let us denote the phase of the Pfaffian by 77. When the configuration is d(> 3)- 
dimensional, we find that 

d n 77 1 d n 

Imlndet7vf = (3-15) 



dAl\ dA% ■ ■ ■ dAH 2 dA%\ dA% ■ ■ ■ dA% 

for n = 1, ■ ■ ■ , (9 — d). This is because, up to (9 — d)-th order of perturbations, 
the configuration stays within 9d configuration and therefore det M remains to be 
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real positive. This means that the phase of the Pfaffian becomes more and more 
stationary as the configuration becomes lower dimensional until one reaches d = 3. 

3.4. Monte Carlo studies of the Euclidean model 

The phase of the Pfaffian makes Monte Carlo studies difficult because of the sign 
problem. However, a new method termed the factorization method 26 )' 27 ) is expected 
to give a definite conclusion on the SSB of SO(10) in the Euclidean model. 

Let us begin by defining the distribution functions for the eigenvalues {Xk} as 



p(xi, ■■■ , xio) = ( n s ( xk ~ Afe v ' ( 3-16 ) 

k 

p^(x lt --- ,x 10 ) = (l[S(x k -X k )) (3-17) 



for the full model and the phase-quenched model, respectively. By definition (3-6), 
these functions vanish unless x\ > ■ ■ ■ > x±q. Applying the reweighting formula to 
the right-hand side of (3-16), one finds that it factorizes as 

p(xi,--- ,xi ) = — p (0) (xi, ■ ■ ■ ,xi ) w(xi, ■ ■ ■ ,xi ) • (3-18) 

The real parameter C is a normalization constant given by*) 

C d =<e* r >o = <cosr>o, (349) 

which need not be calculated in the present method. The function w(xi, • • • , x±o) is 
defined by 

w(xi,--- ,x 10 ) = (e ir ) xij ...,xi = (cosr) xlj ... jX10 , (3-20) 
where ( • ) X1) ..., X10 denotes a VEV with respect to the partition function 

Z X1 ,.., X10 = f dAe~ s " Y[S(x k - A fc ) . (3-21) 

k 

At large N, the expectation values are obtained by maximizing the distri- 
bution function p(x±, • • • , xio) with respect to x\, • • • , x\q. This leads to the coupled 
equations 

d d 

logp (0) (a;i,--- ,xi ) = -t. — \ogw(x 1 , ■ ■ ■ ,x 10 ) for n = 1, • • • , 10 . (3-22) 



The function on the left-hand side and w(xi, • • • , xio) defined by (3-20) can be ob- 
tained by simulating the constrained model (3-21). 

An important observation now 26 ) is that log w (xi, ■ ■ ■ ,x±o) ~ 0(N 2 ) as nat- 
urally expected from the number of degrees of freedom. On the other hand, we 
speculate 78 ^ that the left-hand side (3-22) is O(N) ii]Jx k = (t 2 ) w and x n > r 2 , 



*^ In the second equality, we have used the fact that the phase F flips its sign under the parity 
transformation A\o M> — Aio, which is a symmetry of the phase-quenched model. 
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where I and r are some scale determined by the phase-quenched model. Thus the 
maximum is essentially given by x\ = • • • = Xd = R 2 and Xd+i = • • • = x\q = r 2 with 
R d r w ~ d = £ 10 , where 3 < d < 9. In order to determine which d gives the absolute 
maximum of p(a?i,--- , xio), one can neglect the effect of p(°\xi,--- , xio) since it 
is suppressed by 1/N. Due to the property in the previous section, it is expected 
that logw(x\, • ■ ■ ,xio) becomes maximum for the solution with d = 3. Therefore, 
p(xi, ■ ■ ■ ,x\o) becomes maximum also for d = 3. It is anticipated that the results 
of the Gaussian expansion method can be reproduced by Monte Carlo simulation in 
this way. 

3.5. Definition of the Lorentzian model 

In this section we discuss the Lorentzian matrix model as a nonperturbative 
formulation of type IIB superstring theory in (9+l)-dimensional space-time. 19 ) 
Our starting point is the action 20 ) S = Sb + Si , where 

S h = -^tv([A ll ,A,}[A^A^ , 

Si = tr (^(Cn Q/3 [A„^]) , (3-23) 

with A^ (/i = 0, • • • ,9) and W a {a. = 1, • • • , 16) being N x N traceless Hermitian 
matrices. The Lorentz indices fi and v are raised and lowered using the metric 
r\ = diag(— 1, 1, • • • ,1). The 16 x 16 matrices are ten-dimensional gamma matrices 
after the Weyl projection, and the unitary matrix C is the charge conjugation matrix. 
The action has manifest SO(9,l) symmetry, where A^ and & a transform as a vector 
and a Majorana-Weyl spinor, respectively. The Euclidean model, which has SO(10) 
symmetry, can be obtained from this action by the Wick rotation Aq = iA\o, T° = 
—iTiQ. A crucial difference is that the bosonic part of the action in the Euclidean 
model is positive definite, whereas in the Lorentzian model it is 

tr (i^FH = -2tr (F 0i ) 2 + tr (F tj ) 2 , (3-24) 

where F^ v = —i[A^,A u ] are Hermitian matrices, and hence the two terms in (3-24) 
have opposite signs. 

We define the partition function of the Lorentzian model by 

Z = J dAd>Fe iS = J dAe iSh ¥iM{A) , (3-25) 

where the Pfaffian Pf7W(^4) appears from integrating out the fermionic matrices 
Note that in the Euclidean model, the Pfaffian is complex in general, and its phase 
plays a crucial role in the SSB of SO(10) symmetry. 25 ) -27 )' 71 )' 72 ) On the other hand, 
the Pfaffian in the Lorentzian model is real. Therefore, the mechanism of SSB that 
was identified in the Euclidean model is absent in the Lorentzian model. 

In the definition (3-25), we have replaced the "Boltzmann weight" e~ s used in 
the Euclidean model by e lS . This is theoretically motivated from the connection 
to the worldsheet theory. 20 ) The partition function (3-25) can also be obtained 
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formally from pure M = 1 supersymmetric Yang-Mills theory in (9 + 1) dimensions 
by dimensional reduction. Note, however, that the expression (3-25) is ill-defined and 
requires appropriate regularization in order to make any sense out of it. It turns out 
that the integration over is divergent, and we need to introduce two constraints 

^tr(A ) 2 <^tr(^) 2 , (3-26) 
ltr(^) 2 <£ 2 . (3-27) 

This is in striking contrast to the Euclidean model, in which the partition function 
is shown to be finite without any regularization. 21 )> 22 ) 

Note that e iSh in the partition function (3-25) is a phase factor just as in the 
path-integral formulation of quantum field theories in Minkowski space. However, 
we can circumvent the sign problem by integrating out the scale factor of A^, which 
essentially replaces the phase e lSh by the constraint S*b ~ 0. (Such a constraint is 
analogous to the one that appeared in the model inspired by space-time uncertainty 
principle. 79 )) Without loss of generality, we set L = 1 in (3-27), and thus we arrive 
at the model 

Z = J ^5 (Itr PfA^ ^ , 

(3-28) 

where 9(x) is the Heaviside step function. Since the Pfaffian Y'iM.(A) is real in the 
present Lorentzian case, the model (3-28) can be studied by Monte Carlo simulation 
without the sign problem.*) Note that this is usually not the case for quantum field 
theories in Minkowski space. 

3.6. Monte Carlo studies of the Lorentzian model 

We perform Monte Carlo simulation of the model (3-28) by using the RHMC 
algorithm 38 ^ as in the Id SYM case discussed in section 2. 

In order to extract the "time evolution", we diagonalize Ao, and define the 
eigenvectors \t a ) corresponding to the eigenvalues t a of Aq (a = 1, ■■■ ,JV) with 
the specific order t\ < ■■■ < tjy. The spatial matrix in this basis (t a \Ai\tb) is not 
diagonal, but it turns out that the off-diagonal elements decrease rapidly as one goes 
away from a diagonal element. This motivates us to define n x n matrices A\ ab \t) = 
{t v+a \Ai\t v+b ) with 1 < a, b < n and t = ± J2a=i W« for i/ = 0, • • • , (N - n). These 
matrices represent the 9d space structure at fixed time t. (This point of view can 
be justified in the large- N limit, in which more and more eigenvalues of Aq appear 
around some value t within a fixed interval 5t.) The block size n should be large 
enough to include non-negligible off-diagonal elements. In Fig. 10 (Left) we plot 
the extent of space R(t) 2 = ^tr^4j(i) 2 for N = 16 and n = 4. Since the result 
is symmetric under the time reflection t — > — t as a consequence of the symmetry 

*-* Strictly speaking, the Pfaffian can flip its sign, but we find that the configurations with 
positive Pfaffian dominates as N is increased. Hence, we just take the absolute value of the Pfaffian 
in actual simulation. 
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Fig. 10. (Left) The extent of space R{tf with N =16 and n = 4 is plotted as a function of t for 
five values of k. The peak at t = starts to grow at some critical k. (Right) The 9 eigenvalues 
of Tij(t) with N — 16 and n = 4 are plotted as a function of t for k = 4.0. After the critical time 
t c , 3 eigenvalues become larger, suggesting that the SO (9) symmetry is spontaneously broken 
down to SO (3). 



Aq —t- — Aq, we only show the results for t < 0. There is a critical re, beyond which 
the peak at t = starts to grow. 

Next we study the spontaneous breaking of the SO (9) symmetry. As an order 
parameter, we define the 9x9 (positive definite) real symmetric tensor 

Tij(t) = ^tT{Ai{t)Aj{t)} , (3-29) 

which is an analog of (3-5) in the Euclidean model. The 9 eigenvalues of Tij(t) are 
plotted against t in Fig. 10 (Right) for re = 4.0. We find that 3 largest eigenvalues of 
Tij(t) start to grow at the critical time t c , which suggests that the SO (9) symmetry 
is spontaneously broken down to SO(3) after t c . Note that R(t) 2 is given by the sum 
of 9 eigenvalues of Tij (t) . 

It turned out that one can remove the infrared cutoffs re and L in the large- iV 
limit in such a way that R(t) scales. This can be done in two steps, (i) First we 
send re to oo with N as re = f3 N p (p ~ \)- 80 ^ The scaling behavior is clearly seen 
in Fig. 11 (Left). The scaling curve of R(t) one obtains in this way depends on 
0. (ii) Next we send (3 to oo with L. The two limits correspond to the continuum 
limit and the infinite volume limit, respectively, in quantum field theory. Thus the 
two constraints (3-26), (3-27) can be removed in the large- iV limit, and the resulting 
theory has no parameter other than one scale parameter. 

Let us discuss the second limit (ii) in more detail. We find that the inequality 
(3-27) is actually saturated for the dominant configurations. Therefore, one only 
has to make the rescaling A^ i— >■ LA^ in order to translate the configurations in the 
model (3-28) as those in the original partition function. It turns out that R(t) for 
the rescaled configurations scales in (3 by tuning L and shifting t appropriately. In 
order to see this, it is convenient to choose L so that R(t) at the critical time t = t c 
becomes unity, and to shift t so that the critical time comes to the origin. Then 
R(t) with increasing f5 extends in t in such a way that the results at smaller \t\ scale. 
This is demonstrated in Fig. 11 (Right), where we find a reasonable scaling behavior 
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Fig. 11. (Left) The extent of space R(t) 2 for k = /3N 1/4 is plotted for N = 8, 12, 16 with = 2. 
We plot the results against the shifted time t — t c in units of the size of the universe R(t c ) at 
the critical time. (Right) Similar plot for fixed N = 16 with k, — 2,4, 8. 



for N = 16 with k = 2.0,4.0,8.0. Note, in particular, that the extent of "time" 
increases as k is increased, which is not the case in the bosonic model. 80 ) Thus, 
supersymmetry is crucial for the "emergent time" in the Lorentzian matrix model. 

3.7. The mechanism of SSB in the Lorentzian model 

The SSB of SO(9) looks mysterious at first sight, but we can actually understand 
the mechanism quite intuitively. Let us consider the case in which k is large. Then 
the first term of (3-24) becomes a large negative value, and therefore the second term 
has to become large in order to make (3-24) zero as required in (3-28). Due to the 
constraint -^tr (Ai) 2 = 1, however, it is more efficient to maximize the second term 
of (3-24) at some fixed time. The system actually chooses the middle point t = 0, 
where the suppression on Ai from the first term of (3-24) becomes the least. This 
explains why the peak of R(t) at t = grows as we increase k. 

Let us then consider a simplified question: what is the configuration of Ai which 
gives the maximum -^tr (F^) 2 with fixed -^tr (Ai) 2 = 1. Using the Lagrange multi- 
plier A, we maximize the function G = tr(Fjj) 2 — Atr(ylj) 2 . Taking the derivative 
with respect to Ai, we obtain 2 [Aj, [Aj,Ai\] — XAi = 0. This equation can be solved 
if Ai = x^i for i < d, and Ai = for d < i < 9, where Lj are the representation 
matrices of a compact semi-simple Lie algebra with d generators. Clearly d should 
be less than or equal to 9. It turns out that the maximum of -^tr (Fij) 2 is achieved 
for the SU(2) algebra, which has d = 3, with Li being the direct sum of the spin-^ 
representation and (N — 2) copies of the trivial representation. This implies the 
SSB of SO(9) down to SO(3). The SSB can thus be understood as a classical ef- 
fect in the k — > oo limit. When we tune k with increasing N as described above, 
quantum effects become important. We have confirmed 80 ) that the n x n matrix 
Q = Yli=i A~i(t) 2 has quite a continuous eigenvalue distribution, which implies that 
the space is not like a two-dimensional sphere as one might suspect from the classical 
picture. 
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§4. Summary 

We have discussed the origin of space-time from the viewpoint that matrices 
are the fundamental degrees of freedom in superstring theory. We have seen that 
matrices indeed provide an appropriate description of the quantum space-time at 
the singularities that appear inside a black hole or at the beginning of the Universe. 

In particular, Monte Carlo studies of U(iV) SYM have deepened our under- 
standing on the gauge/gravity duality considerably. In the DO-brane case, it would 
be interesting to investigate whether the duality holds including the 1/N correc- 
tions, which should be compared with the string loop corrections on the gravity side. 
The D3-brane case, which corresponds to the 4d N = 4 SYM, is challenging, but 
we hope that more nontrivial tests of the AdS/CFT correspondence are possible by 
measuring non-circular Wilson loops 17 -* and four-point functions. 81 ) As a new direc- 
tion, Monte Carlo simulation of matrix models obtained by applying the localization 
technique 58 ) to supersymmetric gauge theories is expected to be useful in extending 
the first-principle studies of the gauge/gravity duality to many more examples. 82 ) 

The results for the Lorentzian matrix model, on the other hand, suggest that 
(3+l)-dimensional expanding universe emerges dynamically from type IIB super- 
string theory if the theory is treated nonperturbatively. This may be contrasted 
with the quantum cosmology in the early 80s that aimed at describing the birth of 
the universe 83 ^ within the mini-superspace approximation.*) Note also that the pic- 
ture suggested here is quite different from that in (perturbative) superstring theory, 
where space-time with various dimensions can be obtained by compactification or 
by using D-brane backgrounds. 

The rapid expansion of the three-dimensional space observed in Monte Carlo 
simulation may be interpreted as the beginning of inflation. It would be interesting to 
investigate the microscopic origin of the inflation along this line. Since the mechanism 
of the SSB relies crucially on the noncommutativity of space, it is important to see 
how a commutative space-time appears at later times. Furthermore, it might be 
possible to understand the origin of dark energy found in the present cosmological 
observations and to predict the fate of our Universe by studying the late time physics 
of the Lorentzian model. 

Since superstring theory is not only a theory of quantum gravity but also a 
theory of all the matters and the fundamental interactions among them, it would be 
interesting to see how the Standard Model appears at later times in the Lorentzian 
matrix model. Finding solutions to the classical equation of motions 30 )' 86 )~ 88 ) and 
performing perturbative expansion around them would be an important direction as 
an approach complementary to Monte Carlo simulation. 

We hope that all the ideas and technologies that we learned in Monte Carlo 
studies of QCD would be useful in developing the Monte Carlo studies of superstring 
theory further. 



*' More recently, a nonperturbative approach to quantum gravity has been pursued using the 
causal dynamical triangulation. 31 ' For earlier works that put forward the idea to use matrices for 
cosmology, see Refs. 84). See also Refs. 85) for related works on emergent gravity. 
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